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A supersymmetric analysis is presented for the d-dimensional Dirac equation with central poten- 
tials under spin-symmetric (S(r) — V(r)) and pseudo-spin-symmetric (S(r) = — V(r)) regimes. We 
construct the explicit shift operators that are required to factorize the Dirac Hamiltonian with the 
Kratzer potential. Exact solutions are provided for both the Coulomb and Kratzer potentials. 
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I. INTRODUCTION 

Exact solutions of the Dirac equation that are allowed by certain soluble potentials have always been of interest 
in relativistic quantum theory 1-9]. More particularly, certain aspects of deformed nuclei have been studied for over 
thirty years [ToL lllj ] by means of spin-symmetric and pseudo-spin-symmetric concepts. Ginocchio [l2l - [l3 | showed 
that spin symmetry occurs when the difference between the vector potential V(r) and scalar potential S(r) is a 
constant, V(r) — S(r) = const.; and pseudo-spin symmetry occurs when the sum of the vector potential V(r) and 
scalar potential S(r) is a constant, V(r) + S(r) = const. In the spin-symmetric limit the Dirac equation corresponds 
7-H ■ to a Schrodinger equation which possesses SU(3) symmetry [121 ]. The pseudo-spin symmetry constraint implies a 
degeneracy of the single-nucleon doublets: this can be shown explicitly in terms of the non-relativistic quantum 
numbers (n, £, j = £ + 1/2) and (n— l,£ + 2,j = £ + 3/2), where n, £ and j are the single-nucleon radial, orbital, and 
' total angular-momentum quantum numbers, respectively. 

In this paper, we formulate the spin-symmetric and pseudo-spin-symmetric problems generated by the Dirac equa- 
tion in such a way that they are amenable to an analysis in terms of factorization and shape-invariance methods [15| ] . 
In particular, we provide a supersymmetric analysis of such Dirac Hamiltonians with Coulomb and Kratzer potentials. 
This leads to expressions for all the bound states and corresponding energy eigenvalues for these spin-symmetric and 
pseudo-spin-symmetric problems. In particular we show that the spin and pseudo-spin symmetry limits lead respec- 
tively to distinctly different features for the ground state, which, to our knowledge, have not been noted before. Some 
earlier exact solutions of the Dirac equation within the context of supersymmetric quantum mechanics may be found 
K> ! in Refs. [IM|. 

We consider a single particle that is bound by attractive central vector and scalar potentials, respectively V and S, 
in d > 1 spatial dimensions and obeys the Dirac equation. For central potentials in d dimensions the Dirac equation 
can be written [19( in natural units h = c = 1 as 

i— = HV, where H = > a s p s + (m + S)/3 + V, (1) 
ot ' 

s—l 

m is the mass of the particle, V(r) and S(r), r = |r|, are the spherically symmetric vector and scalar potentials, 
and {a s } and /3 are Dirac matrices, which satisfy anti-commutation relations; the identity matrix is implied after 
the vector potential V. For stationary states, algebraic calculations in a suitable basis lead to a pair of first-order 
linear differential equations in two radial functions {?/>i(r), ip2(f)}- For d > 1, these functions vanish at r — 0, and, 
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for bound states, they may be normalized by the relation 

oo 

(V>i,Vl) + W>2,V>2)= (^l{r)+^l{r))dr = l. (2) 



We use inner products without the radial measure factor n ' because the factor r 2 is already built in to each 
radial function. Thus the radial functions vanish at r — and satisfy the coupled equations 

Ei> x = (V + m + S)iJi + (-d + kd/r)i>2 (3) 
E^j 2 = {d + k d /r)^ + (V -m-S)iP 2 , (4) 

where k\ = 0, kj = r(j + ^j 2 -), d > 1, and r = ±1 and 9 represents the operator 9/cV. We note that the variable r 
is sometimes written ui, as, for example in the book by Messiah [20| . The quantum number fc^ is related to I and J 
for the spin-symmetric and pseudo-spin-symmetric cases as follows: 

fe = {-(^),;-«"- 

and 



h *=u\h: h 7=1-1 (5) 



^ = \ {l 7~A\ lz7r (6) 
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Here i — I + 1 is called the pseudo orbital angular momentum [3 ■ The radial functions are often written tpi = G 
and ij)2 — F, as in the book by Greiner [2l|. We shall assume that the potentials V and S are such that there 
are some discrete eigenvalues E^ dn and that Eqs.Q are the eigenequations for the corresponding radial eigenstates. 
Here n = 0, 1, 2, . . . enumerates the radial wave functions for a given kd- In this paper we shall present the problem 
explicitly for the cases d > 1. Similar arguments go through for the case d — 1: in this case k\ = 0, the states can be 
classified as even or odd, and the normalization ((2} becomes instead (tfjf (x) + ip% (x)) dx = 1. 



II. SHAPE INVARIANCE AND SUSY PARTNER HAMILTONIANS 

A supersymmetric (SUSY) Hamiltonian satisfies the following graded Lie algebra [22j 



where the supercharges Q and commute with H and are nilpotent operators: (Q^) 2 = = Q 2 . Here, the operators 
are given by 

A± = ±±+W(r), (7) 
where W(r) is the superpotential. The associated SUSY partner Hamiltonians H\ and H 2 have the standard forms 

H 1 = ~+V 1 (r), H 2 = -^+V 2 (r). (8) 
dr z dr 

In the case of unbroken SUSY, there is a remarkable aspect of SUSY-QM, namely, except for the zero-energy eigenstate, 
the SUSY partner Hamiltonians Hi and H 2 are found to be exactly isospectral. Here Vi(r) and V 2 (r) are corresponding 
partner potentials satisfying 

V 1 (r;a 1 ) = V 2 {r;a 2 ) + R{a 1 ), a 2 = /(01). (9) 

where ai and a 2 are constants and the remainder R(ai) is independent of r. Often the relation between the potentials 
is given by V 2 {r\a{) = Vi(r;a 2 ) + R(a%) in the literature, but we prefer to use the notation of Sukumar [24|, which 
allows ipi and ip 2 , respectively, to be retained as the upper and lower radial functions in the Dirac spinor. By using 
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the shape invariance condition (|9|) . the entire spectrum of H2 can be found. Thus one can construct the sequence of 
Hamiltonians by using iteration [15j . to find 



j2 s- 1 



H s = --^ + V 2 (r;a s ) + ^R(a k ), a. = f- l {ai), a = 1,2,..., (10) 

k=l 

where f s {a) means s repeated applications /(/(. . . (/(a)) . . .) of the function /. In this fashion, the entire spectrum 
of eigenenergies for the initial Hamiltonian H2 can be obtained algebraically by 

n 

e^ = o, 4 2) = ]>>K). (11) 

fe=i 

Here the superscript (2) denotes the eigenvalue for the Hamiltonian H2- The corresponding eigenfunctions for the Hi 

and H2 are written simply i/ji and %p 2 , and also, more fully, for the corresponding n th excited states, ipn and i/4 . 

(2) 

The unnormalized energy eigenfunction ip„ for the Hamiltonian H2 reads 

^ ex {A+A+...At_ l )^\r-a n + 1) (12) 

but one can also use [l5| 

^\r;a 1 )=A+(r;a 1 )^l 1 (r;a 2 ). (13) 

The complete spectrum is given by 

E%=E£). (14) 
The reader may wish to look in Ref. [l5[ for more details. 

III. FACTORIZATION AND SHAPE INVARIANCE FOR THE COULOMB PROBLEM 

A. Spin-symmetric problems 
In this case ([1]) can be written in the form: 

"(t)-U- E m t E )(t) M> 

in which 



H=^-l + U, (16) 



dr 

1 is the identity matrix, and the matrix U is given by 

/ h& 
' = {-2V(r) -4? 

We shall now use the Coulomb potential V(r) = — ^ in this matrix equation. We can diagonalize the matrix U by 
means of the similarity transformation D~ 1 UD, where 



D 



fed 

V 

1 1 



We define tpi and ^2 as the transformed radial components obtained by 
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If we multiply (|15|) by D 1 from the left, we have 



A+i> 1 = (m + E)^j> 2 (17) 

ltd 



and 



A^^ 2 =[{m + E)^--^{m-E))i> 1 , (18) 



where the operators Aq have the form: 



v 



AZ=±- + f - im + E)-. (19) 



Thus we may now write the eigenvalue equations for ipi and ip2, namely 

Aq A+$i = efa (20) 
AqAq^ = eip 2 , 

where the eigenvalue e is given by 

t 
k 2 

Using (|17l) and (ITBl . we can construct the partner Hamiltonians Hi and H2 

H 1 = A A+, H 2 = A+A^ (22) 



[m + E)% -(m'-E?)). (21) 



+ _£l_± d _ 2 (m + E)- + (m + E) jo 

dr 2 r 2 r hi 



= -^2+ r2 - -2(m + E)- + (TO + ^r^ (23) 
ff 2 = + - 2(to + 25)H + (m + ^) 2 S" ( 24 ) 



By comparing Eqs. (|2"3")l and (j2"4l with Eq. ©, it is clear that ai and a 2 are given by ai = kd and a 2 = k d + 1. By 
using (fTT]l. (fl7]l and ((lSj), we find 



_ m 2 + _ ( ^ n + m)2 = (m + En) 2 y 2 |_ _ ___ j f (25) 

and the entire spectrum can be obtained as 

1 - 1,2 

i4 2 >=£„ = m- <*ffi_ n = 0,l,2,... (26) 

Eq. (f2"U| shows that AJAq and have the same spectrum except when Aq -0 2 = 0. Thus, the ground-state wave 

function becomes 

= r ^ e -(™+ E )ti r (27) 

and the ground- state energy is 

,(2) 1 fe 7 



^=m— — f. (28) 



1 ■ s 



Thus, the first order differential equation (jTTJ) implies 



ipQ — r e, ~ d / dz z 2kd e *5 . (29) 
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First excited state of rp2 is given by oc AJ^q 2 '. The complete set of solutions can be obtained by using (|13|) and 
v 7 m 2 — E 2 = (m + E)v/{kd + n) from the energy relation (|25l) : 

^) oc e-^ r £*"->■ U^f^) (30) 

where C b n {x) are Laguerre polynomials. 

In order to select square integrable wave functions, we have to distinguish different cases. A solution of the 
differential equation A^-02 = yields (j2"T)l with the ground state energy (|2"5)) . We see that ^q 2 ' satisfies the boundary 
conditions when kd > 0, v > 0. However, if we now look at (|29|) . we see that the n = solution for does not satisfy 
the boundary condition when kd > 0, v > 0. Thus, we see that E = Eq is not an eigenvalue of ifi because Hi has 
no corresponding L 2 solution. This means that when spin symmetry occurs, it is not possible for n = to construct 
a normalisable spinor in terms of tp\ and 1^2 • O n the other hand, in the Dirac equation for a central Coulomb field 
problem, without the restriction of spin symmetry, one can indeed find a singlet state when n = 0, for example if 
S = Q. 



B. Pseudo-spin-symmetric problems 



In this case, the matrix equation becomes 



:0 + (t j r£ ) )(£)-U-*"r)(£i- 

Again, the matrix which diagonalizes the matrix with terms = in pip is given by 

D = ( 1 ^ 
^0 1 

If I? -1 is applied to (|3~Tj) from the left, as above, then we obtain 

= (m + E - Hy(m - £) j ^2 (32) 

and 

A-ip 2 = (E-m)iP 1 , (33) 

where the operators can be defined as 

.4- d kd , , u , . 



For this case, the eigenvalue equations are 



where 



Now the partner Hamiltonians become 



A-A+^ = e$ x (35) 



e= (i; 2 -m 2 + ^2(£;-TO) 2 ) . (36) 



d 2 fc d (fc d + l) u 2 u 2 

ff l = + - m )~ + - m ) fc| ( 3? ) 

d 2 fcd(fcd ~ 1) ,« , ,„ ,2^ /ocA 

# 2 = + 2(£-m)~ + (£-m) (38) 
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The parameters are again a\ = k d and a 2 = kd + 1. By using (1361) . we can obtain the spectrum for this case as 

4 2) = -m ]~^ . (39) 

From d32]) and ([33]), we have V>° = °> thus 

^(2) = r fc, e -(B-m)^r j ^(1) =0 (4Q) 

can be obtained. If we look at ip ' in (|^D|) . it is clear that the positive values for kd and negative values for v must 

be used in order to obtain normalisable solutions "tp^ ■ F° r n — 0, v < we have normalisable solutions ■ 

~ (2) 

And the solutions ip\ are given by 

$2) r ^-^f^ c ^-x ( 2 ^^r) . (41) 

IV. THE KRATZER POTENTIAL 

We shall use the Coulomb problem of the previous section as a guide in our analysis of the Kratzer case. We write 
the Kratzer potential]!!! in the form: 

V(r) = ±-~ + c. (42) 

A. Spin symmetric case, S — V 

In the spin symmetric case, we can give the eigenvalue equation by using ([3]) and ((4]): 
d 2 k d (k d + 1) + 2\(E + m) 



( ., , 2(E + m)- + 2{E + m)c) Vi = (E 2 - m 2 )^. (43) 



Now, we introduce the operators as 



4=±^ + A(r). (44) 



Using the Coulomb case of section 3 as a guide, we adopt an ansatz for f3%(r) of the form 



/3iM = ---, (45) 
r Si 



where s% and a± are the positive constants. Partner Hamiltonians are given by 

Hi -A A a - --^ + - 2 — + — 



9 

r sf 



(46) 



and 



1 , d 2 Si (s-i — 1) 2ai a? 

^2 = 4Mo = -— + 11 \ > - -1 + -i. (47) 



<ir 2 



The relation between the potentials is given by 



2 2 



^1 (r, Sl ;v)=V 2 (r, Sl + l;v) + ^~ 1 . (48) 

Here, we present the expressions for the first-order operators as 

A+Vi = mV>2, ^0^2 = (49) 
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If we compare (|46 p -(|48 p and (j43| . we obtain s±, ct\, and the product fiis (to be factorised later): 



and 



Using ([TT]) and (gSJl, we have 



si = i + ^(fc d + i) 2 + 2A(£; + m), ai = »(^ + m), (50) 



99/ s v 2 (E + m) 2 , , 

/xi/ = £ 2 - m 2 - 2(E + m)c + — > — g '—. (51) 



= v 2 (E n + mf ( -1 - * ) . (52) 



This leads to the following general energy relation, which is consistent with results obtained by different methods and 
reported in [25j . 



(n + 1/2 + v/(fc d + 1/2) 2 + 2\{E n + m))y/m 2 - E 2 + 2c(E n + m) = v(E n + to). (53) 
Now we factorize the \iv expression in the form 

fi = m + E (54) 

v 2 

v = E - m-2c+ —{E + m). (55) 

The ground state wavefunction then reads 



^ = r Sl e ~ (56) 



and the ground state energy is 



1 v 

— ~z 2c 

E =m ^-H T . (57) 

And the solutions can be obtained as 

fm+slt. _ . / (rn -X- F,\u \ 

r I . (58) 



^)cxr-e-^ £ 2 — f 2 (" + S)« 



The boundary conditions imply v > 0. Again, i/^ 1 -* does not satisfy the boundary condition of square-integrability. 



B. Pseudo-spin symmetric case, S = —V 
For S — — V, we can give the eigenvalue equation by using ([3]) and (@|: 



v 



4- 2 + fcrf(fcrf ' 1)+ 2 2A(g ' m) - 2(E - m)- + 2(E m)c ) ^ = (E 2 - m 2 )^. (59) 
ar r z r 



And we write the ansatz for /^M as; 



p 2{r)= S -l-^l. (60) 
r s 2 



This leads to partner Hamiltonians H\ and Hi given by: 



H x = A;A$=* + »^-*Z + * (61) 
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dr 2 

We again compare (j59")l and (|62l) and obtain: 

1 



H 2 ~A A - -— + — — + -J- 



s 2 = - + V(kd ~ 1/2) 2 + 2X(E - m), a 2 = v(£7-m). (63) 



Following the same steps as in the spin symmetric case, the energy relation is found to be 



v(E n - m) = (n + 1/2 + v 7 ^ - 1/2) 2 + 2A(£„ - m)) V2c(£ n - m) + m 2 - £ 2 (64) 
and the ground state wavefunction is 

^ 2) =r s *e- {E - m) ^ r (65) 

and the same conclusion is valid for (1651) as given in the Coulomb potential case. Here, if v < we have normalisable 
solutions for the ground state with ip^ = and given in (|65p . Thus, solutions are given by 

^ oc r ^ E ~ m ^ r (2 {E ~™ )V r^ . (66) 

V. SHIFT OPERATORS AND THE KRATZER POTENTIAL 

Introducing an operator Oe, 

d 2 

E = -r 2 -r^ + (2c(E + m) - E 2 + m 2 )r 2 - 2v(E + m)r (67) 
dr z 

we can express (|43|) in the form 

O e ^ Ekd = ~{k d {k d + l) + 2X(E + m))^ Ekd . (68) 
In order to factorize (|6"5|) . let us first use a — ^ E +m)v anc ^ 2c — E + m = Using these parameters, we have 

/ E + m \' 

n = v{ - . 69 

We observe that, (|6"5|) can be obtained from ([2"5| . For the bound states for this potential, we introduce an operator 

O n M 

d 2 2 r 2 

°n = - r2 T^-- r + —> (7 °) 
dr a an 

O n + R ± = Q? l±1 Qt, (71) 

where R± are constants. Our aim is to construct the in order to factorize O n . We assume an ansatz for Q„, 
namely 

Q^ = ±nz^---+n, z = -, [z,p z ] = i. (72) 
dz a n 

If we insert (|72"|) in (jTTj) . we obtain R± = n(n ± 1). However, by use of (|72l QniiQn leads to 

O n + n(n ± 1) = (tt^- - — + n ± l) [±r^- - — + n J . (73) 
\ dr an / \ dr an J 

This means that we have to find an operator such that [26[ 



^ = r— — D ± p P = — ivD* n^l (74) 
nil n 
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and 

(Z>±)t£>± = l, (D^ = DT ±1 . (75) 
With the aid of D^ ±1 D^ = 1, one obtains 

Q* ±1 Qt=O n +n(n±l), n?l. (76) 
Here, by using Q„ can be written as 

Qt = D±T±, T± = ±rj- - — + n. (77) 

ctr an 

From ()76[) . we conclude that |ni?>=e ± |n± 1,£ >. By using the normalization condition, we get 

\e ± f=<n£\(T 7 tyT±\nl>. (78) 

From (|H7J), (EH) and [Hi, we see that 

(^=7^1, <r>=^(3n 2 -^+l)), < rp r >= I, (79) 

and in turn we obtain 



e ± = 



(n± 1) 

>-(n(n± !)-£(£ + !)). (80) 



It is clear that the minimum value for n is n m i n — £ + 1 from (|5D|) and 1 1 >= 0. Thus, if we apply Q^ +1 Q^ +2 ---\^ + 

1 , £ > , we get 

n = N + £+l, N = 0,1,... (81) 

If we use (|8Tj) , (|69|) and ^(fc^ + 1) + 2X(E + m) = £{£ +1), we can obtain the energy relation of the Kratzer potential 
(|53[) for the spin symmetric case. Similar steps can, of course, be used for the corresponding pseudo-spin-symmetric 
problem. 

VI. CONCLUSION 

In this paper we have used shape- invariant techniques to effect a supersymmetric analysis of the Dirac equation with 
the Coulomb and Kratzer potentials, for both the spin-symmetric (S — V) and pseudo-spin-symmetric (S = —V) 
cases. It is shown that the state ip^' is missing for the spin symmetric case; this does not happen with pseudo- 
spin symmetry. Thus the spin-symmetric and pseudo-spin-symmetric cases are qualitatively different. For bound 
states, with either of these potentials, we find that the Coulomb coupling v must have the appropriate sign: for spin 
symmetry, v > 0; for pseudo-spin symmetry, v < 0. This is consistent with the results of Ref. [25|. We have shown 
that the factorization of the Dirac equation for the Kratzer potential can be performed by following similar steps to 
those used for the Coulomb case. To illustrate the significance of this important alternative algebraic approach, we 
have employed the raising and lowering shift operators to obtain the spectrum generated by the Kratzer potential for 
the spin-symmetric case. 
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